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Abstract 

The dynamic structure factor of a normal Fermi gas is investigated by using the moment method 
for the Boltzmann equation. We determine the spectral function at finite temperatures over the 
full range of crossover from the collisionless regime to the hydrodynamic regime. We find that the 
Brillouin peak in the dynamic structure factor exhibits a smooth crossover from zero to first sound 
as functions of temperature and interaction strength. The dynamic structure factor obtained using 
the moment method also exhibits a definite Rayleigh peak (w ~ 0), which is a characteristic of the 
hydrodynamic regime. We compare the dynamic structure factor obtained by the moment method 
with that obtained from the hydrodynamic equations. 
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I. INTRODUCTION 



When discussing collective modes in quantum many-body systems at finite temperatures, 
there are two regimes of interest: the collisionless (or mean- field) regime and the coUisional 
(or hydrodynamic) regime. The mechanisms for the occurrence of collective modes in these 
two regimes differ critically. One of the collective modes is the first sound in the coUi- 
sional regime, which is due to local equilibrium. In a normal Fermi system, zero sound is a 
characteristic sound mode in the collisionless regime. It propagates due to the mean-field 
interaction. This zero sound was first predicted by Landau [1] based on Fermi liquid the- 
ory {2!, and it has been studied in many fields of physics including low-temperature physics, 
nuclear physics (hot nuclear matter), and astrophysics (neutron stars). 

Collective modes with time- dependent density disturbances have been investigated. The 
crossover between first sound and zero sound was first observed by Keen et al. in mea- 



surements of the acoustic impedance between liquid He and a quartz crystal [3|. Density 
fluctuations have been excited in ultracold atoms by deforming the trapping potential j^, 5|. 
An alternative way to probe collective excitations is to use scattering of lights or particles. 
This involves measuring the spectral function of the system, such as the dynamic structure 
factor or the density response function. 

In neutron scattering experiments, the roton spectrum of liquid '^He was observed through 
the spectral function. In early experiments on liquid '^He, several problems were encoun- 
tered in neutron scattering measurements of the dynamic structure factor including a high 
neutron absorption cross-section and interaction between spins. Nowadays, these problems 
have been overcome and the neutron scattering experiments on liquid ^He have also been 
performed . The dynamic structure factor of ultracold atoms has been studied by two- 
photon Bragg spectroscopy in a condensed Bose-Einstein gas and the Bogoliubov spectrum 
was obtained This technique has also been used recently for a Fermi gas [9]. 

In 1958, Abrikosov and Khalatnikov conducted the first theoretical study of the dynamic 
structure factor of a normal Fermi liquid [l^ . They proposed using light scattering to observe 
the zero sound in liquid ^He. The spectral function has also been studied in connection 
with evaluating the Landau parameter by sum rules 
of hot nuclear matter have been studied by taking only two moments: the density and 
current [l^ . In the field of the ultracold quantum gases, the spectral function of normal 
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Fermi gases has been studied 

Urn- 

The dynamic structure factor of normal Fermi gases in both the colhsionless and hydro- 
dynamic regimes has been discussed in detail 16|. By employing the random- phase approx- 
imation, the dynamic structure factor can be discussed even beyond the phonon regime in 
the colhsionless regime. However, the hydrodynamic regime cannot be investigated using 
the same theoretical framework. In this regime, the hydrodynamic equations, which are not 
valid in the colhsionless regime, can be used to calculate the density response function. 

Although the crossover between zero and first sound modes has been extensively studied 
for a long time, there has been no comprehensive study of the dynamic structure factor over 
the full crossover range from the colhsionless regime to the coUisional regime within a single 
theoretical framework. Furthermore, the dynamic structure factor for the crossover from 
zero to first sound has not been explicitly calculated and it is not obvious how it varies as a 
function of temperature and interaction strength. It is thus important to study the dynamic 
structure factor at finite temperatures from the colhsionless regime to the hydrodynamic 
regime within a single theoretical framework. 

In the present paper, we study the dynamic structure factor of a normal Fermi gas over 
the full crossover range from the colhsionless regime to the hydrodynamic regime within 
a single theoretical framework, namely the moment method. The moment method can be 
used to perform systematic analysis and it yields important physical insights. Guery-Odelin 
et al. first applied the moment method to study the collective mode in a trapped classical 
gas [l7|. We recently used the moment method to study excitation spectra of a normal Fermi 
gas in a uniform system ISj; the results for both frequency and damping of the collective 
mode clearly show the crossover from the zero to first sound mode. In the present paper, 
we extend the study of Ref. [l^ by investigating the dynamic structure factor over the full 
crossover range from zero to first sound. 

This paper is organized as follows. Section |Tl] presents the moment method for a normal 
two-component Fermi gas. Section UTTl examines the dynamic structure factor obtained using 
the moment method and discusses the crossover from the colhsionless to the hydrodynamic 
regime. Section IIVI compares the spectral function obtained by the moment method with 
that obtained using the hydrodynamic equations. We focus on the Brillouin peak in the 
dynamic structure factor, which is associated with the sound mode, and the Rayleigh peak, 
which is associated with the thermal diffusion mode. We discuss the results in Sec. El 
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Section |Vl] presents the conclusions. 



II. MOMENT EQUATION AND DYNAMIC STRUCTURE FACTOR 

We start with the following Boltzmann equation: 

dfa{p,r,t) p 

TT. + — ■ Vr/a P, r, t) 

at m 

- Vrf/.(r,t) ■ Vp/.(p,r,t) =Xcoii[/.], (1) 

where the subscript a = {t,]-} represents the spin component. We consider a normal Fermi 
gas with two spin components in the symmetric configuration A^^ = N^. We also assume 
that atoms with different spins collide with an s-wave scattering length of a. The effective 
potential Ucr{r,t) is the sum of the mean-field interaction gn^fj{Y,t) and the external field 
?7ext(i', where no-(r, t) is the local density and g is the interaction strength, which is given 
hj g = Anh'^a/m. In this study, we consider a spin-independent external field. 

We linearize the distribution function around the static equilibrium (denoted by /°(p, r)) 
using /o.(p,r, t) = /°(p,r) + 6f^{p,r,t). It is convenient to write fiuctuations in the distri- 
bution function around the static equilibrium in terms of the average additional energy of 
the particles i/o-(p, r, t), which is defined by 6fa{p,r,t) = {df^/de'^)i'„{p,r,t). We apply a 
relaxation time approximation to the collision integral Icon[fa]- Using this approximation, 
the collision integral can be reduced to 

Xcoii[/.] = -l^^ = --^Su^, (2) 

where r is the relaxation time, fa is the distribution function in local equilibrium, and SUf^ is 
given by 6^^ = u„ — [Aa + B ■ p + Cp^]. z^o-.iocai = A^ + 'B-p + Cp^ is the solution for the local 
equilibrium. The coefficients Ao-, B, and C are determined below using the conservation 



18| as the relaxation time r in Eq. 



law. We use the viscous relaxation time given in Ref. 

We now consider an external field with the form [/ext(r,t) = f/cxt(q, i^)e*'''''"~'^*^ This 
leads to a plane-wave solution of the linearized Boltzmann equation, which is represented 
as z/o-(p, r, t) = z/o-(p, q, a;)e*'^'' '"~'^*\ and 5na-{r,t) = 6na-{q,,u))e^^'^''^^'^*\ The linearized Boltz- 
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mann equation for z/o-(p) is now given by 
5/° 



p • q 



m 



P ■ q. 



1^ 



[z/.(p)- (A^ + B-p + C/)] 



(3) 



where we omit q and uj in v^j and (5no- for simplicity. 



We expand the fluctuation in terms of spherical harmonics as Vaiv) = 

oo I ^ 

S ^ z/™;(p)P;™'(cos 6')e*™'^. Multiplying Eq. (E]) by e"*™ ''^ and integrating it over 0, we 

l=Qm=-l 

obtain the reduced form of the linearized Boltzmann equation: 



1=0 



dfa (pq 



95° 



cos 6') [gSn^^ + f/ext(q, w)]^-, 



m,0 



1^ 



z/r(p)^r(cos 0) - (A, + B • p + Cj9^) <5™,o 



«=0 



(4) 



We see that only the mode m = is coupled to the external potential, and thus we take the 
mode m = which corresponds to the longitudinal wave. 

When we take moments corresponding to the number of particles, the momentum, and 
the energy, the collision integral vanishes due to the conservation law. The coefficients Ag-, 
B, and C in the relaxation time approximation are determined from these conservation laws. 
The resultant equations are: (z/o-,o) — A^Wafi — CWa,2 = 0, ^ {{p^(t,i) — BWa,2) = 0, and 
Yl {{P^^ufl) - A^W^,2 - CW^^i) = 0. The function W^^n is defined as 



dp ^ndf. 



(27rn)3^ del' 

and we assume longitudinal sound (i.e., B || k). We defined the moment {p"'i^a,i) as 

dp df^ 



(5) 



{P^^^,i) 



(27rn)3 del 

where we used the notation i^^j"^{p) = i'a,n{p)- 

Multiplying Eq. by p^Pi{cos6), where P;"^=°(cos 



(6) 



Pi(cos6'), and integrating it 
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over p, we obtain the moment equation, which we summarize as: 



2/ - Im 



(7) 



2L + 6 m m 



I 

T 

il 

re 



Wo 



1 

e 



wjvl_ 

Wo' 

Wn+2W2 



Wo 

W2Wn 



Wn+2W2 

Wo 



T 2W2 



Wo 



iy„+i(5z,i?7ext(q,c<;), 



m 



(8) 



W. 



where we define 6 = 2 {W4, - W2 /Wo). We used W„ = W„^n = yy-a,n, assummg a 
population-balanced gas. In the absence of an external field (f/ext = 0), solutions of Eq. 
(jH]) give the frequency (or sound velocity) and damping of normal modes. Detailed calcu- 



18|. From the solution of Eq. 



lations of the normal-mode solutions are presented in Ref. 
([8]) including f/ext (q,w), we can calculate the density 6ntot = {^0), which can be written in 
terms of the density response function as (5ntot(q, 1^) = i^)f4xt(q, i^)- 

The density response function yfq, u) is related to the dynamic structure factor ^(q, w) 
by the dissipation fiuctuation theorem [19]. Using the detailed balance condition ^(q, w) = 
e^'^S{q, —u), and the relation Imx(q, w) = — 7r[S'(q, w) — ^(q, —oj)], the dynamic structure 
factor can be written as 



^(q,^) 



1 



71 1 



lmx(q,a;). 



(9) 



We now discuss collective modes that use the dynamic structure factor ^(q, u) calculated 
from the moment equation. The present moment method re pro duces the excitation only in 
the phonon regime q < kp (as discussed in a previous study [l8|), and hence our results are 
valid for phonon regimes. 



6 



III. DYNAMIC STRUCTURE FACTOR FROM THE COLLISIONLESS TO HY- 
DRODYNAMIC REGIME 

A. Strong coupling case 

Figure [1] shows the spectral function in a strongly coupled system as a function of uj 
and T. We used the renormalized frequency and temperature uj/{vpq) and T/Tp, where 
vf = {h/m)(37c'^Ntot/VY^^ is the Fermi velocity and Tp = E-p/k-Q is the Fermi temperature. 
Ey is the Fermi energy, which is given by E^ = mVp/2, A'tot is the total number of particles, 
which is given by A^tot = + N^, and V is the volume. The dynamic structure factor is 
normalized as follows: S{ijJ,T) = S{ijj,T)VeF/Ntof We choose the parameters q = O.OS/cp 
and a = 15, where a = gNtot/iVep) [19], and take moments up to / = n = 30. Figured] 
also shows the sound velocity calculated by the moment method (solid line), the zero sound 
(dashed line), and the first sound (dotted line) in the u-T plane. 

It shows that there is a smooth crossover between zero and first sound in the structure 
factor. As in Ref. [l8|, the sound eigenmode exhibits a crossover between zero and first 
sound. Correspondingly, the peak in the structure factor associated with the sound mode 
(the Brillouin peak) transitions smoothly from the zero sound mode to the first sound mode 
with increasing temperature. As expected, the Brillouin peaks are quite sharp in both the 
collisionless and collisional regimes, where the damping rate is small. The peak in the 
crossover regime from zero to first sound is quite broad, reflecting a short lifetime. The peak 
width in the hydrodynamic mode is discussed in Sec. llVi 

Figure [2] shows the spectral function as a function of frequency u and coupling constant a. 
We set the parameters as g = O.OS/cp and T = 0.025Tf and take moments up to / = n = 31. 
The peaks in the collisionless and collisional regimes are quite sharp. The Brillouin peak 
becomes broader in the crossover regime from zero to first sound. The zero and first sound 
modes have very similar velocities in the strong-coupling regime, which makes it difficult 
to distinguish between the two regimes solely based on the mode frequency. However, the 
reduction in the intensity of the dynamic structure factor with increasing coupling strength 
a is a clear indicator of a crossover between zero and first sound. 
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FIG. 1: The dynamic structure factor S{uj, T) as a function of uj and T, calculated by the moment 
method. In the oj-T plane, the solid line denotes the phase velocity calculated by the moment 
method [l^], the dashed line gives the phase velocity of zero sound, and the dotted line gives that 



of first sound. 

B. Weak-coupling case 

The situation is more complicated for a weakly coupled system. It is difficult to distin- 
guish collective modes from other excitations when the moment method is used to determine 
the eigenvalues of a normal Fermi system because the eigenvalue for coherent oscillation 



overlaps with those for the particle-hole continuum. Because of this complication, Ref. [18 1 
omitted detailed discussion of the collective mode in a weakly coupled system with changing 
temperature. We now study collective excitations in a weakly coupled system in terms of 
the dynamic structure factor. Figure [3] shows a plot of the spectral function of a weakly 
coupled system. We choose the parameters as g = O.Olfcp and a = 1 and take moments 
up to Z = 77, = 31. For reference, the zero sound (dashed line) and first sound (dotted line) 
frequencies are plotted in the u-T plane. 

The collective modes in the coUisionless and coUisional regimes are easily found from 
the spectral intensity peak. The first sound mode clearly has a single broad peak at high 
temperature. Zero sound can be also identified as it has a narrow, high-intensity peak 
slightly above the Fermi velocity at low temperature. However, in the crossover regime from 
zero to first sound, the intensity is low and it is affected by other modes. 
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FIG. 2: The dynamic structure factor S{uj, a) as a function of uj and a, calculated by the moment 
method. The sound velocity calculated by the moment method (solid line) [18] is in the w-a plane. 
The sound velocities of zero sound (dashed line) and first sound (dotted line) are also plotted. 

Thus, it is difficult to distinguish the collective mode in the crossover regime from zero 
to ffist sound from the structure factor in a weakly interacting system; this aspect differs 
from the strong coupling case. Detailed behavior of the spectral intensity in the crossover 
and coUisionless regimes below uo < qvp ( the region corresponding to the particle-hole 
continuum) is affected by the number of moments that we take for numerical calculations. 
However, the main results noted above remain the same. 



C. Rayleigh peak 



In classical hydrodynamics, the response function also includes a peak due to the thermal 
diffusion mode at u 







16 



20 



211 ]. which known as the Rayleigh peak 



expression for the hydrodynamic response function is given in Refs. 



An explicit 



20, \2M j it is discussed in 

the next section in the context of the moment method. As discussed in Ref . [l^ , the normal 
mode solutions in our moment method include the thermal diffusion mode. The dynamic 
structure factor calculated for Eq. (jlj) involves the associated Rayleigh peak. Figure H] shows 
a plot of the spectral function over a wide range in the u-T plane using the same parameters 
as those in Fig. [H The Rayleigh peak appears at w = in the hydrodynamic regime, and 
it appears with increasing temperature. We compare this result with the hydrodynamic 
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FIG. 3: The dynamic structure factor S{uj, T) calculated by the moment method as a function 
of uj and T for the weakly coupling case a = 1. The dashed line shows the phase velocity of zero 
sound and the dotted line shows that of first sound in the oj-T plane. 




FIG. 4: The dynamic structure factor S{uj, T) calculated by the moment method as a function of 
a; and T using the same parameters as those for Fig. [1] The solid line in the oj~T plane represents 
the phase velocities of the sound. There are two peaks: the Brillouin peak at high velocity and the 
Rayleigh peak near a velocity of zero. 

response function in the next section. 
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IV. DYNAMIC STRUCTURE FACTOR IN HYDRODYNAMIC REGIME 



In this section, we briefly review the density response function in the hydrodynamic 
regime and we compare this result with the dynamic structure factor obtained using the 
moment equation. The average additional energy of particles in the local equilibrium i^a-,iocai 
is given by z/o-,iocai = Ao- + B ■ q + Cp^, as mentioned in Sec. [Tll Neglecting the departure 
from local equilibrium 5v^ = v„ — z^o-,iocai on the left-hand side of Eq. ([3]) and also f/ext, we 
obtain 



IT 



p ■ q 



pq 



[A. + B ■ p + Cp"] 



(10) 



Substituting z/g- = '^o-.iocai + ^'^o- in Eq. ([3]) and taking the zeroth, first, and second moments 
of the linearized Boltzmann equation, we obtain the hydrodynamic equations given by 



B q 

3m 



(11) 



> {uj^W^^2-A^ ^q-C ^q 



m 



q 



m 



B ■ q 



3m 



(Vrt,4 + W^i,4) - AiKm'CTq' = 0, 



(12) 



(13) 



where the collision integral is zero in the hydrodynamic regime because of the conservation 
law. K is the thermal conductivity k = —kBf^riWQ — /W9) /(6m^), and r] is the shear 
viscosity 77 = — 2rW4/(15m^). Detailed derivations are given in Ref. |l8 |. 

The above equations can be written in terms of fluctuations in the total density 5ntot = 
EAA<rW^,o + CW,,2], velocity v = -B, and energy 6E = J2JA,W^,2 + CW,,2]/{2m). 
The resultant hydrodynamic equations can be expressed in matrix form: 



K 



q ■ V 



( 



\ 





m 



V / 



(14) 
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where the 3x3 matrix K is given by 



K = 



_g<f_ 

2m 



2W2 
3m 

(jj — i- 



Wa 



Wo 



\ 2mWo 3m2 







(15) 



The parameter 7 is defined by 



7 = 



WoiW^-gWj) 



and the rate is given by 



gWo) 

" " ~ (1^4 - gWi){W4Wo - wl 
Solving the matrix equation ( IT4|) for 5ntot, we obtain 

1 



tot 



detir 



32 



'12-"-33j — <-^cxt 

m 



2W2 2w + ^r«7 



ext • 



(16) 

(17) 

(18) 
(19) 



3m^ detfT 

We neglect the second-order terms in the transport coefficients k and rj since they are small 
in the hydrodynamic regime. This allows the determinant of the matrix K to be reduced to 
detir = {u^ - n^){uj + iV^) + 2iTuj^, where 



mo 



m 



cq 



and 



rjq 



nTq^rn^ 



(20) 



(21) 



W2 {Wi - gWi) 

Taking this determinant to be zero, we obtain the eigenmodes of the hydrodynamic modes 

(to first order in k, and rj): u = ±Q — iT and u = —iTf^ 18|. fl is the frequency of the sound 

mode and T is its damping rate, c in Eq. (l20ll is the sound velocity. The damping rate 

is that of the heat diffusion mode. 

Using the above approximations for the zeros of detK, we can explicitly derive the density 

disturbance induced by the external field in terms of the frequency and damping rates of 

the eigenmodes. This is given as 5ntot(q, w) = x(q, i^)f4xt(q, 1^), where the resultant density 

response function x(q, u) is given by 

2W2 2 uj + iT^j 
2^ 



x(q,w) 



3m^ 



{u - n + iT){u + n + iT){u + iT^) 



(22) 
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FIG. 5: The dynamic structure factor S{u},T) obtained from the hydrodynamic equations given 
in Eq. (I23p . using the same parameters as those for Fig. [TJ 



The absorptive susceptibility can be reduce to 



Imx(q,w) 



2W2 r u;(7-i)r. ^ 

Sm^c^ _|_ ^2 ^ 

(a;2 _ (]2)2 + (2wr)2 



(u;2 _ ^2)2 + (2wr)2 



(23) 



where we used ^ F^, and also i7 ^ F. This absorptive susceptibility has two peaks: the 
Rayleigh peak at a; = and the Brillouin peak at u; = f2. 

Figure [5] plots the dynamic structure factor obtained from the hydrodynamic equations as 
a function of u and T. We used the same parameter set as that in Fig. [TJ Comparing Fig. [5] 
with Fig. m shows that the structure factor calculated using the hydrodynamic equations 
differs from that obtained by the moment method in the crossover and collisionless regimes, 
as expected. The spectral intensity decreases in those regimes because of the high damping 
rates of hydrodynamic modes. In contrast, Fig. H] clearly shows that the moment method 
can correctly describe the dynamic structure factor in both the crossover and collisionless 
regimes, whereas the hydrodynamic equations cannot. 

To evaluate how effective the moment method is in the collisional regime, we perform a 
more quantitative comparison between the results obtained using the moment method and 
those obtained using the hydrodynamic equations. Figure [6] shows a plot of Rayleigh peaks 
for several temperatures as a function of u. The solid, dashed, dotted and dot-dashed lines 
represent results for T/Tp = 0.2, 0.1608, 0.1216, and 8.24 x 10~^, respectively. The thick 
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FIG. 6: Rayleigh peaks for several temperatures as a function of oj. Solid, dashed, dotted, and 
dot-dashed lines represent results for T/Tp = 0.2, 0.1608, 0.1216, and 8.24 x 10~^, respectively. 

lines are the results obtained using the moment method and the thin lines are those obtained 
using the hydrodynamic equations. It is not possible to distinguish them since they exhibit 
identical behavior in the hydrodynamic regime. At lower temperatures, the Rayleigh peak 
becomes too weak to distinguish the two results. 

Figure [7] shows a plot of Brillouin peaks for several temperatures as a function of u. 
The solid, dashed, and dotted lines represent the results for T/Tp = 5 x 10^^, 4.02 x 10^^, 
and 3.04 x 10~^, respectively. The thick lines show the results obtained using the moment 
method and the thin lines show the results obtained using the hydrodynamic equations. As 
the temperature decreases and the system approaches the crossover regime, the difference 
between the dynamic structure factor obtained by the moment method and that obtained 
from the hydrodynamic equations increases. In both Figs. [6] and U\ we used the same 
parameters as in Fig. [1] We conclude that our moment method can reproduce the dynamic 
structure factor in the hydrodynamic regime reasonably well. 

We comment on the Landau-Placzek ratio in our formulation. Comparing Eq. f l23|) with 
the formula for the absorptive susceptibility in terms of the thermodynamic quantities in 



where Cp and Cy are the specific heats at constant pressure and volume, respectively. We 
used ntot = 2no- = —2W2/ {3m) and the relation mc^ = {dp/dntot)\s = {Cp/Cv)idp/dntot)\T, 
where p is the pressure and 5* is the entropy. The Landau-Placzek ratio ^lp, defined as the 



Eq. (4.44a) in Ref. |2oj (also Eq. (87a) in Ref. |2l|), we find that Cp/C^ is equal to 7, 
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FIG. 7: Brillouin peaks for several temperatures as a function of lo. Solid, dashed, and dotted 
lines represent results for T/Tp = 5 x 10"^, 4.02 x 10"^, and 3.04 x 10^^, respectively. Thick lines 
are the results obtained by the moment method and the thin lines are the results obtained from 
the hydrodynamic equations. 



ratio of half the area under the Rayleigh peak to that under one Brillouin peak 20| , is given 
by £:lp = (Cp/Cv) — 1 [22]; hence the Landau-Placzek ratio can be obtained as 



Slp =1 -I 

WoWi - Wi 



(24) 
(25) 



Wiil-gWoY 

At T = 0, we have Xn,a = — 3A^totPF/(4V^£:F), and hence the Landau-Placzek ratio ^lp 
becomes zero. Figure [8] shows that the temperature dependence of the Landau-Placzek 
ratio for a = 15. As mentioned above, this ratio approaches zero at T = and it is a 
monotonically increasing function of temperature. To derive the Landau-Placzek ratio, it 
is usually necessary to calculate the specific heats at constant pressure and volume Cp and 
Cy. However, in our formalism it can simply be obtained from the function Wn- 



V. DISCUSSION 



In experiments of ultracold atomic gases, two-photon Bragg scattering is a well-developed 
tool for studying the dynamic structure factor M. This type of study is analogous to 
the classic work by Abrikosov and Khalatnikov [lO|, who proposed using light scattering 
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FIG. 8: The Landau-Placzek ratio elp obtained from our formalism as a fmiction of temperature. 



to observe zero sound in liquid ^He. With the recent dramatic progress in experimental 
techniques for ultracold atomic gases, the coupling strength a can also be controlled because 
of the Feshbach resonance. The results in the present paper should be observable in principle. 

Hu et al. recently discussed the density response function in superfluid gases in the two- 
fluid hydrodynamic regime |23j] . They showed that the first and second sound have different 
relative weights in the dynamic structure factor S{q, u) and the density response function 
Imx(q, u) 23|. This is clearly seen in the relation between both functions shown in Eq. (|9]). 
When uj/ik-oT) <^ 1, ^(q, ~ — A;B^Imx(q, w)/(7r/jiX'). Hu et al. pointed out that the extra 
factor l/u leads to a large enhancement of the peak associated with the low-frequency second 
sound in ^(q, w) [23|. This is also true for the Rayleigh peak in the present study. Figure [9] 
plots the imaginary part of the density response function — Imx(w, T) as a function of u and 
T using the same parameters as for Fig. [H The density response function is normalized by 
xioj^T) = xiy>,T)V£Y/Niox.- Due to the absence of the factor l/cj, the Rayleigh peak near 
w = cannot be observed. Thus, the Rayleigh peak can be observed through the dynamic 
structure factor, whereas it has a negligibly small weight in the density response function. 
A localized potential that turns off after a short duration, which is used in experimental 
studies of ultracold atomic gases, will also be useful for studyin g th e Rayleigh peak because 
the excitation of density pulses is proportional to Imx(q, c<;)/c(; |23l. l24|. 

The moment method can be applied to other systems, such as transverse zero sound and a 
collective mode in a Fermi gas with dipole interaction. The spectral function of a polarized 
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FIG. 9: The imaginary part of the density response function — Imx(a;, T) as a function of uj and T 
using the same parameters as for Fig. [TJ The solid hne in the uj~T plane shows the phase velocity 
of the sound obtained by the moment method. Only the Brillouin peak is visible at high velocity. 
The Rayleigh peak near zero velocity is not visible. 



normal Fermi gas at the unitarity limit has recently been studied at T = [15|. The 
method described in the present paper can also be used to determine the density response 
of a polarized normal Fermi gas with a strong interaction at finite temperatures. We intend 
to apply the moment method to these systems in the future. 

In the present paper, we applied the moment method to the Boltzmann equation with a 
simple mean-field interaction and calculated the dynamic structure factor in the crossover 
regime between the coUisionless regime and the hydrodynamic regime. Our method pro- 
vides a controlled calculational tool as long as the perturbative control of the interaction is 
possible. In the case of a normal Fermi gas near the unitarity limit, the Boltzmann equation 
(JT]) will be modified by replacing the mean-field interaction Ucr = gn^cr (with g = Anh'^a/'m) 



with the real part of self-energy 25| associated with the many-body T matrix T which is 



calculated through the ladder approximation |26l . 12 7| . For the collision integral in the Boltz- 
mann equation, the scattering cross section da/dfl is also given in terms of the many-body 
T matrix T through da/dn = m'^\T\^ / [A-Kh^f Q-Q- 

The above approach involves difficulty when one solves this Baltzmann equation directly; 
but, one can still apply the moment method if one uses the so-called unitarized vacuum 
scattering matrix, i.e. the T matrix neglecting the effects of the medium given hYT = g/ (1 + 
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iap/h), where p is the relative momentum of the scattering particles 27H31|. This approach 
using the vacuum scattering is valid for the high temperature regime (T — Tc)/Tc > 1 291. 



The effects of the medium becomes significant at lower temperatures {T — Tc)/Tc < 1 281.129 1 



We note that the 
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3^. As 



because of the phase transition associated with the Cooper instability 
regime T < 2Tc near the unitarity limit corresponds to the pseudogap regime 
a result, the moment method will be applicable for the normal Fermi gas near the unitarity 
limit at higher temperature than that for the pseudogap regime. Outside the perturbative 
regime, the results may be trusted qualitatively but only to the same extent that one can 
trust qualitative descriptions of the mean-field theory for static equilibrium properties. 



VI. CONCLUSION 



We studied the spectral function of a normal Fermi system at finite temperatures from 
the coUisionless to the hydrodynamic regime. We solved the Boltzmann equation accurately 
using the moment method and we determined the dynamic structure factor in the crossover 
and coUisionless regimes as well as in the hydrodynamic regime as a function of the tempera- 
ture and coupling strength. We compared the results obtained by the moment method with 
those obtained using the hydrodynamic equations in terms of the Rayleigh and Brillouin 
peaks. We also briefly commented on the Landau-Placzek ratio. 

It has been generally difficult to study the dynamic structure factor in the crossover 
regime between the coUisionless and hydrodynamic regimes because of its complexity. This 
study describes a powerful method for calculating the dynamic structure factor over the 
whole crossover range, from the coUisionless to the hydrodynamic regime, within a single 
framework. 
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